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1. Introduction. The celebrated idea by Nambu |l| to replace the fundamental Higgs 
boson in the standard model (SM) by a top quark condensate, which is responsible for 
dynamical symmetry breaking, has been investigated in some detail in refs. 0. In frames 
of the Nambu- Jona-Lasinio mechanism 0, and under some conditions, it was shown that in 
the large- A^ limit the physics of such Nambu- Jona-Lasinio (NJL) model is equivalent to the 
physics of the SM with elementary scalar fields. 

However, the original NJL model [0] with gauge fields and some of its modifications are 
considered now as non-renormalizable effective theories where the presence of an ultraviolet 
cut-off A at the loop diagrams is a necessary condition. It is clear that there are different 
possibilities to extend the NJL model. Among such possibilities, a quite interesting one is 
the introduction in the original Lagrangian of higher-derivative terms in the four fermion 
interaction [Q |^ or in the kinetic term |^. It may be shown |^ that the physics of such 
generalized NJL models is still equivalent to the physics of the SM. Moreover, the additional 
advantage appears, that the ultraviolet cut-off A shows up at the Lagrangian level. That may 
help discuss the ambiguities which turn out in the choice of a cut-off 0, ^ Q and violation 
of background gauge invariance. From another side, the inclusion of higher derivatives in 
effective theories gives the possibility to take into account the structural effects of the medium 
and external fields, like external gauge and gravitational fields or condensates. 

It is the purpose of this work to consider a higher derivative four-fermion model in 
curved spacetime and to study the influence of an external gravitational field to the effective 
potential of composite scalars in the large--^ expansion. In particular, we discuss the phase 
structure and the dynamically generated fermion mass at non-zero curvature. This can be 
important in early Universe considerations. Notice that similar questions for the standard 
NJL model in curved spacetime have been already investigated in refs. [0, |^, ^ (for a study 
of X(p^ theory see |TT]]). In particular, it was shown that curvature may induce first order 



phase transitions |]^, |^. We will here discuss whether these considerations go over to the 
higher derivative four fermion model at issue. 

2. The model and its effective potential. We start by presenting the model which we 

set out to study (it was introduced in [Q). Its Lagrangian density in curved spacetime (for 
an introduction to quantum field theory in curved spacetime, see [|^ is given by 



C = ^P^J'^{x)V,iJ + -^^ \ Ai (^^)' + 3A2 



^(1-2^)^ 



(1) 



where Nc is the number of fermionic species, 'y^{x) the Dirac matrices in curved space, A 
a cut-off parameter, and Ai and A2 are coupling constants. We work in the -^^-expansion 
scheme. 



By introducing some auxiliary fields xi ^^^ X2 we can give a description of this non- 
renormalizable theory by means of the action 



S 



d^x^\^i-f''{x) 



v.ij - iVeA' [y^+^]- xii^i^ + V3x2^ (1 - 2^] V 



• (2) 



Our purpose is to study the influence of external gravity on the dynamical breaking and 
restoration patterns of the symmetry possessed by the Lagrangian (^, which is given by 
the transformations "0 —>■ j^ip and i/j —>■ "075 . If we refer to the action (0), the symmetry 
is realized by adding the following transformations for the auxiliary fields: xi ~^ ~Xi ci-nd 

X2 -^ -X2- 

The effective action of this model in the Nc -^ 00 limit is given by 

Ai A2, 
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e// 



AM^ + ?l+^le// 



(3) 



with 



i {Voiy^ Tr \nliy7-m + 2^3x2^ J 



where m = Xi + v^X2 • Thus, to leading order in the ^^-expansion we have the gap equations, 
as follows 

dVleff -1^ 1 
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(4) 
(5) 



First of all, we have to calculate the fermionic propagator, which satisfies 

li](/-m + 2v^X2^ J Qix, x') = Six, x') . 

Once this propagator is obtained (apart from terms which disappear under the Tr operation) 
it will be immediate to produce equations (^) and (|^) in explicit form and also the dependence 
of the effective potential itself on the fields Xi ^"^^ X2- The details of this rather lengthy 
derivation are written in the Appendix, where we give details about the calculation of the 
effective potential up to terms linear in the curvature. 

The final outcome is shown below in natural variables, which are obtained by using the 
following definitions 



R 



Xi 



X2 



r = — ■, xi = — , X2 = — , a = xi + ^3x2 , v 
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The 'dimensionless' effective potential v is given by 
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Note that in deriving this expression we have taken into account only terms up to quartic 
order on the fields a and X2- The gap equations — which one obtains by differentiating v 
with respect to xi and X2 and equating the results to zero — are, respectively 
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3. Symmetry breaking in curved spacetime. We can study now the influence of gravity 
on the symmetry breaking pattern of the theory. A simple inspection reveals that a positive 
curvature tends to protect the symmetry of the vacuum, and also a negative one may trigger 
the breaking of the symmetry. 



In flat spacetime the symmetry is broken whenever either Ai or A2 are greater than ^'^ 



2 

A2-! 



that is why we shall take in the sequel k\ and k2 to be l\ — %'n'^ and I2 — 87r2 respectively. 
We may illustrate these remarks by studying the evolution of the minimum, given by two 
coordinates (which we choose to be a and X2) in several circumstances. Before including 
gravity, it is worth discussing in more detail the situation in flat space. In fact, several cases 
may be analyzed (see p). To introduce this study it is convenient to define 



±i,\ ^ k- 



57r' 



A 



where Af means Aj > 87r2 (for the + sign) or Aj < 87r2 (for the — sign). One sticks here 
to the case when \ii <^ A. We only repeat the two situations given when ^ > 3^1 and 
{Xi , X2), or //2 < 3//i and (A^^, A2"). In both cases the results may be summarized by 
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X2 = -V^xi — 
^2 ^ /^?A*i 
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As a first example, consider the case in which the couphng constants are such that the 
symmetry is not broken in fiat space, that is to say, Aj < Stt^. Let us see how the situation 
is modified as we move from negative to positive curvature. In all the figures, we have 
represented the values at the minimum of a and X2-, which correspond to the dynamically 
generated mass (m) and X2 transformed to natural units, as given in expression (|^). We 
observe that there is a negative value of the curvature above which the symmetry is restored 
(see Figs. 1 and 2). We observe a continuous phase transition. Some comments about this 
fact are done below. 

We can see the equivalent plots (Figs. 3 and 4) with the only difference that one of the 
couplings is such that the symmetry is broken at R = 0. One can identify a minimum positive 
value of the curvature for the symmetry to be restored. The character of the transition is 
still continuous. 

This observation may be distressing if one recalls the results of Q, Q , where the N JL model 
was studied under the influence of a gravitational field in different situations. The conclusion 
of those works was always that if the coupling constant is greater than the critical value in fiat 
space, there is a first order phase transition at some positive value of the curvature. In fact 
there is no actual contradiction. The point is that, till now, we have concentrated ourselves 
in cases where the coupling constants are around Stt^. We observe that if we explore regions 
in the space of parameters where A2 is much smaller, and Ai is greater than Stt^ — which 
would correspond to a limit where our theory approaches the Gross- Neveu model — then 
there is a first order phase transition for some positive value of the curvature. 

Let us now illustrate what happens if we keep the curvature constant and vary the values 
of the couplings, to study how the domain which guarantees that the symmetry is not broken 
is modified by the presence of curvature. The situation of negative curvature is exemplified 
in Figs. 5 and 6. We see that the symmetry is broken from a negative value of /ci, in other 
words, the region of parameters of the theory where the symmetric vacuum is stable has 
shrinked. 

To finish this discussion we consider the case of positive curvature. Here one can see 
that the aforementioned domain has been enlarged by the influence of positive curvature 
(see Figs. 7 and 8). 

In summary, we have calculated the effective potential of composite scalars in the four- 



ferinion model (Q) in curved spacetiine, and numerically investigated the phase structure of 
the theory. 

There is an interesting question about this model — namely whether it would be possible 
for some generalization of it to be represented in renormalizable form. As we see, if A2 = 
in (|l]), this is perfectly feasible as long as one takes into account those terms in the field xi 
which transform the theory into a renormalizable Yukawa-type model which describes near 
the critical point the physics of chiral symmetry breaking. As we see, the higher derivative 
term in [Q] acts against such generalization at A2 7^ 0. 
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Appendix: Expansion of the effective potential up to linear curvature terms 

In this appendix we give details of the intermediate steps which are required in the 
lengthy calculation of the effective potential in the linear curvature approximation (for the 
representation of propagators in local momentum expansion and weak curvature, see 
The first step is the inversion of the operator given by 

2" 



z y - m + 273x2^) (iW + m- 2v^X2^ j = - ^^ + U - 2v^x2^ j 



(11) 



where □ stands for the spinor d'Alembertian. Henceforth we shall use in our intermediate 
calculations, the relation 

R-a/3^1/ = — 7 7T idafLgpu " gaudfifj.) , (12) 

n [n — I) 

which is valid in a constant-curvature spacetime. As we are restricting ourselves to the 
first nontrivial term in an expansion in the number of derivatives of the metric tensor, this 
apparent simplification will not introduce errors to this order. Again, keeping only terms 
proportional to the curvature, we may write 

R. / ,, ,, o ,,,,\ ^ ^ R ^,1 \ ^ 2 R 

n {n 

being the y'^ Riemann normal coordinates, and 



a = d'+ /" [y^y^ - y\^^) d,d^ + —^^J\y'd, + \-y'd^ , (13) 

3?2 m — 1 ^ ^ n[n — \] in 
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n (n 
Now we write the equation satisfied by the inverse of the operator given in (0) 

{A,-A2 + ...)G{y) = 5{y) 

The functional G may be given as an expansion in derivatives of the metric tensor G 
Go + G2 + . . ., where 

AGo(l/) = 5(1/), 

G2 = G0.4.2G0 



(15) 



Let us give some exphcit forms of these expressions: 
d" k exp {—iky) 



Go 

^2Go 



d^k 



2' 



(16) 
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Some definitions are in order: 
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Now Q = i3G, where B = i JJ + a — 2-^/3x2^2"' ^^*^ ^^ ^^y expand B = Bq + B2 + ■ ■ ■ 
(again according to the number of derivatives of the metric tensor). The exphcit expression 
of the operators Bq and B2 are: 

Bo = z^ + a-2V3x2^, (18) 



n 



where the symbol " means that the Dirac matrices involved satisfy the anticommutation 
relations {7", 7^} = "if^ . 

We need only to expand the expressions shown below 

^0 = SoGo , (19) 

Q2 = B2G0 + B0G0A2G0. (20) 

In particular, the explicit form of Qo and B2G0 are 
f dJ'k exp{-iky) 

r rl^-k R / 1 2 1 \ 

^=«» = -/(2^^^P<-^-^)2^^-=F(M;r3T)'>l(^-V5;^h2(A7 + j) ^ (22) 

Equation ( |22D is true, apart from terms which have an odd number of 7 matrices — and 
which will not contribute to the effective potential. Finally, hi and h2 are defined as: 

hl(A;2) = n{n+l)g{P) + 2P{2n + 3)g'{k'^)+ (23) 

+4{ey g"ie) - {Ak'ne) + 2nne)) , 

h2(A;2) = ngi^) + 2k''g\k'') . 

Putting all these pieces together we may construct the propagator Q (apart from terms 
with an odd number of 7 matrices) and perform the momentum integrations ensuing from 
the Tr operations in (^) and (|^). As for (^, remember that yj"^ = D + ^ and □ is given in 

(in. 
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Figure captions 



Fig. 1 The coupling constants are kept constant at tlie values given by ki = —2 = k2. One 
sees a continuous phase transition. 

Fig. 2 The coupling constants are kept constant at the values given by ki = —2 = k2. There 
is a continuous phase transition at some negative value of the curvature. 

Fig. 3 The coupling constants are kept constant at the values given by ki = —1, k2 = 0.2. 
Now we see a continuous phase transition. 

Fig. 4 The coupling constants are kept constant at the values given by ki = —1 k2 = 0.2. 

Fig. 5 The curvature is such that r = —0.0008 and k2 is kept constant at k2 = —0.8. 

Fig. 6 The curvature is such that r = —0.0008 and k2 is kept constant at k2 = —0.8. 

Fig. 7 The curvature is such that r = 0.0008 and ^2 is kept constant at k2 = —0.8. Now 
the range of values of the coupling constants which maintain the symmetry of the vacuum 
is enlarged as the curvature gets positive. 

Fig. 8 The curvature is such that r = 0.0008 and k2 is kept constant at ^2 = —0.8. 
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Fig. 4 
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Fig. 6 
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